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Abstract 

We investigate lensing and waveguiding properties of an atomic Bose- 
Einstein condensate for ultraslow pulse generated by electromagnetically in- 
duced transparency method. We show that a significant time delay can be 
controllably introduced between the lensed and guided components of the 
ultraslow pulse. In addition, we present how the number of guided modes 
supported by the condensate and the focal length can be controlled by the 
trap parameters or temperature. 
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1. Introduction 

Quantum interference effects, such as electromagnetically induced trans- 
parency (EIT) ug, can produce considerable changes in the optical prop- 
erties of matter and have been utilized to demonstrate ultraslow light prop- 
agation through an atomic Bose-Einstein condensate (BEC) This has 
promised a variety of new and appealing applications in coherent optical 
information storage as well as in quantum information processing . How- 
ever, the potential of information storage in such systems is shadowed by 
their inherently low data rates. To overcome this challenge, exploitation of 
transverse directions for a multimode optical memory via three dimensional 
waveguiding of slow EIT pulse j3] has been recently suggested for BECs 
Transverse confinement of slow light is also quintessential for various pro- 
posals of high performance intracavity and nanofiber slow light schemes (See 
e.g. Ref. |6| and references therein). Furthermore, temperature dependence 
of group velocity of ultraslow light in a cold gas has been investigated for an 
interacting Bose gases Q. 

A recent experiment, on the other hand, has drawn attention that ul- 
tracold atomic systems with graded index profiles may not necessarily have 
perfect transverse confinement due to simultaneously competing effects of 
lensing and waveguiding jsj. The experiment is based upon a recoil- induced 
resonance (RIR) in the high gain regime, employed for an ultracold atomic 
system as a graded index waveguiding medium. As a result of large core 
radius with high refractive index contrast, and strong dispersion due to RIR, 
radially confined multimode slow light propagation has been realized jsj]. As 
also noted in the experiment, a promising and intriguing regime would have 
few modes where guided nonlinear optical phenomena could happen [sj . 

It has already been shown that the few mode regime of ultraslow waveg- 
uiding can be accessed by taking advantage of the sharp density profile of 
BEC and the strong dispersion provided by the usual EIT [if. Present work 
aims to reconsider this result by taking into account the simultaneous lensing 
component. On one hand, the lensing could be imagined as a disadvantage 
against reliable high capacity quantum memory applications. Our investi- 
gations do aid to comprehend the conditions of efficient transverse confine- 
ment. On the other hand, we argue that because the lensing component is 
also strongly delayed with a time scale that can be observably large rela- 
tive to the waveguiding modes, such spatially resolved slow pulse splitting 
can offer intriguing possibilities for creating and manipulating flying bits of 
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information, especially in the nonlinear regime. Indeed, earlier proposals to 
split ultraslow pulses in some degrees of freedom (typically polarization), face 
many challenges of complicated multi-level schemes, multi-EIT windows, and 
high external fields Quite recently birefringent lensing in atomic gaseous 



media in EIT setting has been discussed [10j|. The proposed splitting of lens- 
ing and guiding modes is both intuitively and technically clear, and easy to 
implement in EIT setting, analogous to the RIR experiment. 

The paper is organized as follows: After describing our model system 
briefly in Sec. El EIT scheme for an interacting BEC is presented in Sec. 
EJ Subsequently we focus on the lensing effect of the ultraslow pulse while 
reviewing already known waveguiding results shortly in Sec. |H Main results 
and their discussion are in Sec. Finally, we conclude in Section 



2. Model System 

We consider oblique incidence of a Gaussian beam pulse onto the cigar 
shaped condensate as depicted in Fig. ([I]). Due to particular shape of the 
condensate two fractions, the one along the long (z) axis of the condensate 
and the parallel the short (r) axis, of the incident Gaussian beam exhibit 
different propagation characteistics. The r-fraction exhibit the lensing effect 
and focused at a focal length (/) while the axial component would be guided 
in a multi mode or single mode formation. 




Figure 1: Lensing and waveguiding effects in a slow Gaussian beam scheme with an ultra- 
cold atomic system. 

The angle of incidence (8) controls the fraction of the probe power con- 
verted either to the lensing mode or to the guiding modes. When both lensing 
and waveguiding simultaneously happen in an ultraslow pulse propagation 
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set up, an intriguing possibility arises. Different density profiles along ra- 
dial and axial directions translates to different time delays of the focused 
and guided modes. Due to significant difference in the optical path lengths 
of guided and focused components, an adjustable relative time delay can be 
generated between these two components. As a result, these two components 
become spatially separated. The fraction of the beam parallel to the short 
axis of the condensate undergoes propagation in a lens-like quadratic index 
medium resulting in a change in the spot size or the beam waist of the output 
beam. We call this as the lensing effect and the corresponding fraction as the 
lensed-fraction. The fraction of the beam propagating along the long axis is 
propagating in weakly-guided regime of graded index medium that can be 
described in terms of LP modes. This is called as guiding effect and the cor- 
responding fraction is called as guided fraction. We use the usual Gaussian 
beam transformation methods under paraxial approximation to estimate the 
focal length. In addition our aim is to estimate time delay between these two 
fractions. The output lensed and guided fractions of the incident beam are 
delayed in time relative to each other. In general temporal splitting depends 
on modal, material and waveguide dispersions. For a simple estimation of 
relative time delay of these components, we consider only the lowest order 
modes in the lensed and the guided fractions, and take the optical paths as 
the effective lengths of the corresponding short and long axes of the conden- 
sate. In this case we ignore the small contributions of modal and waveguide 
dispersions and determine the group velocity, same for both fractions, by as- 
suming a constant peak density of the condensate in the material dispersion 
relation. 

3. EIT scheme for an interacting BEC 

A Bose gas can be taken as condensate part and thermal part at low 
temperature. Following Ref. density profile of BEC can be written by 
p{f) = p c (r) + pth(r), where p c (r) = [(p(T) - V(r))/U }Q(p - V(f)) is the 
density of the condensed atoms and pth is the density of the thermal ideal 
Bose gas. Here Uq = 47rh 2 a s /m; m is atomic mass; a s is the atomic s-wave 
scattering length. 9(.) is the Heaviside step function and T c is the critical 
temperature. The external trapping potential is V(r) = {m/2){uj 2 r 2 + uj 2 z 2 ) 
with trap frequencies u r ,u z for the radial and axial directions, respectively. 
At temperatures below T c , the chemical potential p is evaluated by p(T) = 
Ptf(Nq/N) 2 / 5 , where ptf is the chemical potential obtained under Thomas- 
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Fermi approximation, /j TF = ((hw t ) / '2) (15 Na s /a h ) 2 ^ , with tu t = (u^cj 2 ) 1 / 3 
and = y/h/ (c^w 2 ) 1 / 3 , the average harmonic oscillator length scale. The 
condensate fraction is given by No/N = 1 — x 3 — sC(2)/C(3)x 2 (l — x 3 ) 2//5 , 
with x = T/T c , and £ is the Riemann-Zeta function. The scaling parameter 
s is given by s = ^ T F/k B T c = (l/2)C(3) 1 / 3 (15iV 1 / 6 a s /a h ) 2 / 5 . 

Treating condensate in equilibrium and under Thomas-Fermi approxima- 
tion (TFA) is common in ultraslow light literature and generally a good ap- 
proximation because density of an ultracold atomic medium is slowly chang- 
ing during the weak probe propagation. The propagation is in the order of 
microseconds while atomic dynamics is in millisecond time scales. Due to 
weak probe propa gati on under EIT conditions, most of the atoms remain 



in the lowest state[12[. The validity of TFA further depends on the length 
scale of the harmonic potential. If the length scale is much larger than the 
healing length, TFA with the harmonic potential works fine. The healing 



length of the BEC is defined as £ = [1/ (87rna s )} 1 / 2 [13j where n is the density 
of an atomic Bose-Einstein condensate and it can be taken as n = p(0,0). 
We consider range of parameters in this work within the range of validity 
of TFA. The interaction of atomic BEC with strong probe and the coupling 
pump field may drive atomic BEC out of equilibrium. This non- equilibrium 
discussion is beyond the scope of the present paper. 

We consider, beside the probe pulse, there is a relatively strong coupling 
field interacting with the condensate atoms in a A-type three level scheme 
with Rabi frequency Q c . The upper level is coupled to the each level of the 
lower doublet either by probe or coupling field transitions. Under the weak 
probe condition, susceptibility x f° r the probe transition can be calculated as 
a linear response as most of the atoms remain in the lowest state. Assuming 
local density approximation, neglecting local field, multiple scattering and 
quantum corrections and employing steady state analysis we find the well- 
known EIT susceptibility [H, EH, Xi^ = r i z i fo r either radial (r) or axial (z) 
fraction of the probe pulse. Total EIT susceptibility for BEC in terms of the 
density p can be expressed as Xi — Pi Xi i n the framework of local density 
approximation. Here Xi is the single atom response given by 

_ M 2 z(-zA + T 2 /2) 

Xl e h (r 2 /2 - *A)(r 3 /2 - iA) + nl/ a' 1 } 

where A is the detuning from the resonant probe transition . For the ultra- 
cold atoms and assuming co-propagating laser beams, Doppler shift in the 
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detuning is neglected, « is the dipole matrix element for the probe transi- 
tion. It can also be expressed in terms of resonant wavelength A of the probe 
transition via \i = 3eoh\ 2, y/S7i 2 , where 7 is the radiation decay rate of the 
upper level. Y 2 and T 3 denote the dephasing rates of the atomic coherences 
of the lower doublet. At the probe resonance, imaginary part of x becomes 
negligible and results in turning an optically opaque medium transparent. 



4. Propagation of beam through a quadratic index medium 

4-1. Lensing effect 

We can neglect thermal part of a Bose gas due to the high index contrast 
between the condensate and the thermal gas background so that p = p c . 
We specifically consider a gas of N = 8 x 10 6 23 Na atoms with = O.57, 
7 = 2tt x 10 7 Hz, T 2 = 7 x 10 3 Hz, and fi c = 27. We take u r = 160 Hz and 
oj z = 40 Hz. For these parameters, we evaluate \' = 0-02 and x" — 0.0004 at 
A = O.I7, where x' an d x" are the real and imaginary parts of x, respectively. 
Neglecting x" 1 the refractive index becomes n = y/1 + x'- I n the z direction 



it can be written as 15, 16 



n(z) = n,[l - fcz 2 ] 1 ' 2 , (2) 

where no = (l+px'1/Uo) 1 ^ 2 and quadratic index coefficient is (3 2 = x'i m u 2 / (ZUqUq) 
Thomas- Fermi radius for the axial coordinate is given by Rtf z = \/2/i(T)/mu;^ 
Expanding Eq. (J2]) in Taylor series, the refractive index reduces to 

n{z)*n [l- l -[}y]. (3) 



The refractive index as a function of z is shown in Fig. [21 With such 
a refractive index profile, the atomic medium can act as a thin lens for the 
component of the probe in the radial direction. In the case of a medium with 
a quadratic index variation, we can estimate the focal length by using geo- 
metrical or beam optics 15|, ll6|. In the paraxial geometrical optics regime, 



where the angle made by the beam ray with optic axis is small, the differ- 



ential equation satisfied by the ray height s(r) is d 2 s/dr 2 + (3 2 r = [15|, [16 
The initial ray height and initial slope are r\ 
solutions of the differential equation are s(r) 



9i = ds/dr\i respectively. The 

= riCos((3 z r) + 1/(3 Z sm(j3 z r)9i 
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Figure 2: z dependence of the refractive index of a condensate of 8 x 10 6 23 Na atoms at 
T = 296 nK under off-resonant EIT scheme. The other parameters used are O c = 1.5 x 7, 
A = 0.1 x 7, M = 23 amu, A = 589 nm, 7 = 2?r x 10 7 Hz, T 3 = O.57, T 2 = 7 x 10 3 Hz. 



and 9(z) = ds/dr = —j3 z sin(/3 2 r)rj + cos(f3 z r)8i. Hence, the ray transforma- 
tion matrix which connects the initial and final ray vectors is given by 



Eq. 01]), where r is the propagation length in the medium. 15|, [l6|, [17 



V -/3,sin(/3,r) cos(/3 2 r) )' 1 1 

The same ray transformation matrix given by Eq. (j4j) can be used to 
investigate the effect of the quadratic medium on the light in the beam optics 
regime. Consider a collimated Gaussian beam incident on a focusing ultra 
cold medium represented by My. The q parameter q(r) of the Gaussian beam 
which describes the spot size and the radius of curvature of the beam is given 
by 

-1 — (5) 



q(r) R(r) nirw(r] 



2 ■ 



where w(r) is the spot-size function, and n is the medium refractive index. 
If the incident beam is collimated, R(r) — > 00 , and the initial g-parameter 
qi becomes iz , where the Rayleigh range z is given by z = tiUq/X . Here, 
Wq is the initial beam waist on the BEC and we assumed that the back- 
ground refractive index is nearly 1. After a distance of r inside the BEC, the 
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transformed g-parameter q 2 will given by 



iz cos(/3 r r) + j- sin(/3 2 r) 
2 -iz /3 z sin(/3 2 r) + cos(/3 2 r) ' 

We have two cases to consider for focusing. If the BEC is very thin, the 
focal length will be approximately given by 1/ (3 z L r where L r is the transverse 
width of the BEC. If the BEC length is not negligible, the next collimated 
beam will be formed inside the BEC at the location Lf, given by /3 z Lf = tt/2. 

4-2. Waveguiding effect 

Along the radial direction, similar to previous treatment, refractive index 
profile in the radial direction can be written as 

where n = (1 + ^x'i/Uo) 1 ^ 2 and f3% = Xi™^ / (2U nl) . Thomas-Fermi radius 
is given by R TFr = y / 2/j l (T)/moj%. 




Figure 3: r dependence of the refractive index of the condensate. The parameters are the 
same as Fig. ([2]) 

The refractive index which is given in Eq. (J7J) is plotted as a function of r 
in Figj3j Expanding Eq. (jjj) in Taylor series, the refractive index reduces to 
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n(r) m n (l — 0.5/3^r^). Such an index behavior is analogous to the one of 
a graded index fiber [5j. Slow light propagation through the condensate can 
be described similar to that of the weakly guided regime of the graded index 
fiber, where the optical modes can be given in terms of linearly polarized 
modes (LP modes). The mode profiles are determined by solving the wave 
equation which reduces to the Helmholtz equation Jsj] . 

We use the cylindrical coordinates as the refractive index n(r) is axially 
symmetric. The wave equation for the axial fraction of the probe field is 
[V 2 + k 2 ]E = 0, where k 2 = k 2 + k 2 and V 2 is the Laplacian operator 
in cylindrical coordinate. Here k r is the radial wave number and k z is the 
propagation constant in the z direction. The solution for the wave equation is 
E = ip{r-) cos(/0) exp[i(ut — k z z)\. Here / = 0, 1, 2, 3, ... and is the absolute 
phase. If we put this solution in the wave equation, we get the Helmholtz 
radial equation 

[d 2 /dr 2 + (l/r)d/dr + p 2 {r)]ip{r) = 0, (8) 

in which p 2 {r) = [k^n 2 ^) — k 2 — l 2 /r 2 ) (l5j |. Here ko can be expressed in 
terms of resonant wavelength A of the probe transition ko = 2n/\ = 1.07 x 10 7 
(1/m). 

We use transfer matrix method developed in Ref. |5( in order to solve the 
Helmholtz equation. Eq.flSJ). 

We assume that the atomic cloud can be described by layers of constant 
refractive index, such that, their indices monotonically increase towards the 
center of the cloud. By taking sufficiently large number of thin layers such 
a discrete model can represent the true behavior of the refractive index. 
For constant index, analytical solutions of the Helmholtz equation can be 
found in terms of Bessel functions. Such solutions are then matched at the 
shell boundaries of the layers. Doing this for all the layers, electromagnetic 
boundary conditions provide a recurrence relation for the Bessel function 
coefficients, whose solution yields the wave number k. The mode profiles 
determined by this method are shown in Fig. (jBJ), and Fig. (J2J). 

5. Results and Discussions 

5.1. Lensing Properties 

The condensate can act as a lens or a guiding medium depending on the 
Rayleigh range of the probe relative to the effective length of the condensate 
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For the cigar shaped BEC geometry we consider, the radial fraction of 
the probe is subject to lensing while the axial fraction is guided. Let us 
first examine the focal length of the condensate for the radial fraction of the 
incident probe. For that aim we need to determine the effective radial length 



d 3 rr 2 p(r, z) 



v 



1/2 



(9) 



corresponding to the radial rms width of the density distribution. Radial 
component of the probe field enters the medium and converges at a focus 
distance, determined 

/ = «b (10) 
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Figure 4: Behavior of the focus length with respect to trap frequency for the radial direc- 
tion and scattering length. The solid curve is for a s = 5 nm and the dashed curve is for 
a s = 7 nm. The other parameters arc the same as Fig. @ 

Focal length is plotted in as shown in Fig. 01]) as a function of trap 
frequency for the radial direction and with the scattering length at a tem- 
perature T = 296 nK. The increase in the radial trap frequency leads to 
decrease in the radial size of the BEC, which causes the rise of the focal 
length seen in Fig. (j3J). In contrast, the size of the condensate increases 
by atom-atom interactions, characterized by the scattering length a s , which 
causes an overall reduction of the focal length. 




10 



5.2. Waveguiding Properties 

For the cigar shaped BEC geometry we consider again the axial fraction 
of the probe is subject to guided. The effective axial length is determined by 



47T 

TV 



poo poo 

/ rdr / dzz 2 p(r, z) 
Jo Jo 



1/2 

(11) 



The axial length L z is an effective length corresponding to the axial width 
of the density distribution. Group velocities (v g i, i — r,z) of the different 
density profiles for both radial and axial directions can be calculated from 
the susceptibility using the relations 

1 1 7r<9Yj 

= I ±L (Y)\ 

v gi c A dA ' 1 ; 

Here, c is the speed of light and imaginary part of the susceptibility is neg- 
ligibly small relative to the real part of it. EIT can be used to achieve 
ultraslow light velocities, owing to the steep dispersion of the EIT suscepti- 
bility x [sl- I n general temporal splitting depends on modal, material and 
waveguide dispersions. For a simple estimation of relative time delay of these 
components, we consider only the lowest order modes in the lensed and the 
guided fractions, and take the optical paths as the effective lengths of the 
corresponding short and long axes of the condensate. In this case we ignore 
the small contributions of modal and waveguide dispersions and determine 
the group velocity, same for both fractions, by assuming a constant peak den- 
sity of the condensate in the material dispersion relation. Density profiles 
along radial and axial directions lead to different time delays of the focused 
and guided modes. Due to significant difference in the optical path lengths 
of guided and focused components, a certain time delay can be generated 
between these two components. 

Taking into account transverse confinement, ultraslow Gaussian beam 
propagation along the long axis of a BEC can be described in terms of multi- 
ple LP modes, propagating at different ultraslow speeds \^ . In the geometry 
and system we consider here, for instance, there can be found approximately 
44 modes at temperature T = 42 nK. We compare the delay time of the 
slowest mode (the lowest guided mode) with that of the lensing mode which 
can be respectively calculated by toz = L z /v gz and tpr = L r /v gr . Thermal 
behaviors of the time delays are shown in Fig. ([5]). Fig. ([5]) indicates that 
due to different refractive indices seen by the axial and radial fractions of the 
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Figure 5: Thermal behavior of the time delays of the guided and focused fractions. The 
solid and dashed lines are for toz and tor , respectively. The parameters are the same as 
Fig. ([!]) except O c = O.57 

probe, significant temporal delay, ~ 50 //s, can occur between them at low 
temperatures T ~ 42 nK. Time delays for each fraction decreases with the 
temperature. Besides, the relative difference of their time delays diminishes 
with the temperature. The condensed cloud shrinks due to the increase in 
temperature, so that effective lengths of the BEC diminish. Therefore, just 
below T c , time delays drop to zero. 

At low temperatures waveguiding modes can occur in an atomic Bose- 
Einstein condensate because of an atomic Bose-Einstein condensate. Trans- 
lation of the Thomas- Fermi density profile of the condensate to the refractive 
index makes the mediunigain waveguiding characteristics analogous to those 
of a graded index fiber p, |8] . For the same set of parameters, for which the 
radial fraction of the probe undergoes lensing effect, the axial fraction is 
guided in multiple LP modes, as the corresponding Rayleigh range is larger 
than the effective axial length. The lowest two LP modes are shown in Fig. 
OH]), and Fig. ([7]). The total number of modes that can be supported by 
the condensate is determined by the dimensionless normalized frequency V 
which is defined as V = (u/c)R(nl — l) 1 / 2 where R = y^2jj{T) JnuJ% 
The radius of a condensate that would support only single mode, LP o, is 
found to be R ~ 1 /.mi. In principle this suggest that simultaneous lensing 
and single LP mode guiding should also be possible. 
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Figure 6: Contour plot of intensity of the V'oo exp[i(wi — /3 z z)] (LP o mode). The param- 
eters are the same as Fig. |5) 



6. Conclusion 

We investigate simultaneous lensing and ultraslow waveguiding properties 
of atomic condensate under EIT conditions by using a quadratic-index model 
of the medium. The focus length, relative time delay, and multiple guided 
mode characteristics are determined taking into account three dimensional 
nature of the system. In particular, dependence of focus length on atom- 
atom interactions via s-wave scattering length, and on trap frequency and 
temperature are examined. Our results reveal how to select a suitable set of 
system parameters to tune Rayleigh range and the aspect ratio of the cloud 
to make either lensing or guiding more favorable along a particular direction. 
We have shown that the focus length can be calibrated by the transverse 
trap frequency, incoming wave length, scattering length and temperature. In 
addition, time-delayed splitting of ultraslow Gaussian beam into radial and 
axial fractions is found. 
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